Time-of-flight imaging method to observe signatures of antiferromagnetically ordered 

states of fermionic atoms in an optical lattice 
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We propose a simple method to detect the antiferromagnetic (AF) state of fermionic atoms in 
an optical lattice by combining a time-of-flight (TOF) imaging method and a Feshbach resonance. 
In this scheme, the nontrivial dynamics of fermionic atoms during the imaging process works as a 
probe with respect to the breaking of the translational symmetry in the AF state. Precise numeri- 
cal simulations demonstrate that the characteristic oscillatory dynamics induced by the scattering 
process that transfers an AF ordering vector appears in TOF images, which can be easily observed 
experimentally. 



PACS numbers: 05.30.Fk, 37.10.Jk, 71.10.Fd 



Cold atoms trapped in optical lattices are ideal model 
systems for investigating the many-body problems that 
have been studied for many years in condensed matter 
physics [l|, [3- These systems are highly controllable 
and thus considered to be quantum simulators. In fact, 
the quantum phase transition from a metallic state to a 
Mott insulating state has been successfully realized using 
fermionic ''"K atoms in optical lattices 0,|3|- The antifer- 
romagnetic (AF) transition in the optical lattice system 
is now a major concern for condensed matter physicists. 
The thermodynamic properties of the AF states in opti- 
cal lattices have been studied in detail in the recent the- 
oretical work [3, [a] . On the other hand, no experimental 
method has yet been established for detecting such AF 
states, and this is an important issue. Several proposals 
have been reported recently based on the Bragg diffrac- 
tion of light \J\ , the detection of noise correlations @ , and 
high-resolution in situ imaging !9!], however these meth- 
ods require extremely sophisticated experimental tech- 
niques. This convinces us of the need for complemen- 
tary easy-to-use approaches employing conventional ex- 
perimental methods. 

In this paper, we propose a simple method to detect AF 
states by combining time-of-flight (TOF) imaging and a 
Feshbach resonance. A TOF imaging method is a useful 
standard tool for investigating the optical lattice systems 
[lOl lll| . Using this method. Kohl et al. observed the 
Fermi surface of two-component fermionic atoms in an 
optical lattice, and also observed certain excitations of 
atoms to higher Bloch bands in combination with a Fes- 
hbach resonance [l4|. Our present study is highly moti- 
vated by this report. We carry out the precise numerical 
simulations and then demonstrate that the nontrivial dy- 
namics of fermionic atoms during a TOF imaging process 
associated with a Feshbach resonance is sensitive to the 
breaking of the translational symmetry of the system. 
We discuss the use of this feature as a probe for the AF 
state. 

We consider a virtual system consisting of fermionic 
'"^K atoms with two different hyperfine states in an opti- 
cal lattice. We assume that the trapping potential in the 



vertical direction is much larger than that in the hori- 
zontal direction, which effectively makes the system two- 
dimensional (2D) 13|. The parameters are set as the 
relevant values as in the experiments |12| : trapping po- 
tential frequency Wtrap = 27r x 150 Hz, scattering length 
asc — 1.5 nm, lattice laser wavelength A = 826 nm, lat- 
tice laser beam waist w — SO /im, and lattice potential 
depth Vo = UEr, where Er = h'^/{2M\^) is the recoil 
energy and M is the mass of "^^K atoms. Furthermore, 
we consider the balanced population of each hyperfine 
state. 

The present system can be well described by the 2D 
Hubbard model with a harmonic trapping potential Il4| . 
The Hamiltonian is written as, 



T~i — 2^ 2^ Jacl^a^^r+Aaa + 2_^ Z_^ C^Q^^rta 
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where clg.^{cr(7a) creates (annihilates) a fermionic atom 
of the ath Bloch band with pseudospin (hyperfine state) 
"■(— tii) on the lattice site at position r, and rira-a — 
'^laa'^rcra- Here, Jq is the nearest-neighbor hopping inte- 
gral, Ea is the onsite potential, Vta is the curvature of the 
harmonic trapping potential for atoms in the ath band, 
Uap is the on-site interaction between two atoms in the 
ath and the /3th bands with different pseudospins, and fi 
is the chemical potential. With the exception of n, these 
parameters are uniquely calculated from the setup pa- 
rameters mentioned above, Vq, A, and so on. The chemi- 
cal potential fi is determined by fixing the total number 
of atoms N. The gap between the first and second Bloch 
bands A = £2 ^ £i is very large in our setup, and the 
population of the higher bands is very small. This allows 
us to carry out the present numerical calculation within 
the second Bloch bands. In the following, we describe for 
simplicity the quantities of the first Bloch band without 
indices, for instance, Ji = J and Un = U . 

There are four energy scales in the present system; 
temperature T, the bandwidth W = 8 J, the interaction 



strength U and the trapping energy Et = VtN/2'K Q. 
Using these scales, we can regard the calculated results 
as being independent of the details of the parameters, 
such as Vt, ^, and the system size [3, @]. Here, we set W 
as an energy unit. In Ref. [6|, we investigated the AF 
transition temperature of a 2D optical lattice system. We 
found that the highest transition temperature is obtained 
at around U/W ^ 1.5 and Et/W ^ 0.5. In the following 
calculations, we always realize U/W ^ 1.6 and Et/W ^ 
0.5 by choosing the setup parameters. 

Before explaining the details of our calculations, it 
is useful to summarize how we obtain TOF images of 
fermions in an optical lattice. As mentioned in Refs. jl3 |. 
first the lattice potential is linearly ramped down, namely 
Vo is linearly decreased to zero. This ramping down 
procedure induces nontrivial dynamics, because the pa- 
rameters of the Hamiltonian ([T]) change as Vq decreases. 
Then, the trapping potential is turned off suddenly, and 
accordingly unconfined atoms expand ballistically in sev- 
eral milliseconds. A TOF image of the atoms is finally 
obtained with a CCD camera. Such a TOF image ideally 
corresponds to the momentum distribution of the atoms. 
However, the experimentally observed image must be af- 
fected by the nontrivial many-body dynamics during the 
lattice ramping down. In addition, TOF images are cer- 
tainly infiuenced by the ballistic expansion process in a 
finite time and the resolution of the CCD cameras [15| . 
Our numerical simulations deal with all these effects. 

Here we briefly explain procedures of our calculations. 
First, we calculate the initial thermal states using a con- 
ventional mean-field approach. The nontrivial dynamics 
during the lattice ramping down is studied in the same 
mean-field manner. We decompose the ramping time 



ramp hi^to Small Segments characterized by the short time 



scale St{= iramp/i)- At each segment, we recalculate the 
parameters in the Hamiltonian ([1]) according to the time 
dependence of Vq, and simulate the time evolution dur- 
ing 6t using the exact diagonalization. We confirm that 
L is sufficiently large and that this decomposition does 
not affect the results. Next we simulate the ballistic ex- 
pansion of atoms in the time ^tof using the obtained 
time-evolved states. A TOF image is finally evaluated as 
the spatial distribution of expanded atoms. We assume 
that the CCD is a 120 x 120 array of 2 ^m square pixels so 
that the intensity of a TOF image I{x, y) corresponds to 
the number of expanded atoms in a 2 /im square around 
the position {x, y). 

We show the properties of the initial thermal states. 
In Fig.[IJa) and (b) , we show the number of atoms with 
t pseudospin for the different temperatures T /W ^ 0.8 
and T/W ^ 0.3, respectively. As shown in Fig.[IJa), the 
paramagnetic metallic (PM) state appears at the higher 
temperature T/W ^ 0.8. In Fig.[ljb), we find a checker- 
board pattern around the center of the trapping poten- 
tial, meaning the AF state appears at the lower temper- 
ature T/W ~ 0.3. 

Next, in Fig.[IJc) and (d), we show TOF images cal- 
culated for these two different initial states. Here, we 
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FIG. 1: (Color Online) Number of atoms with t pseudospin 
(nrt) at (a) T/W ~ 0.8 and (b) T/W ~ 0.3. Time-of-flight 
image obtained after 2 ms linear lattice ramping down at (c) 
T/W ~ 0.8 and (d) T/W ~ 0.3. 
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FIG. 2: (Color Online) (a) Schematic diagram of the sequence 
of the TOF imaging to detect the AF state. TOF images 
taken by the present sequence for (b) the PM state (obtained 
at T/W ~ 0.8) and (c) the AF state [T/W ~ 0.3). The image 
size is 240 /im square. 



set t 



ramp 



=2ms and <tof =8ms. After 8 ms of ballis- 
tic expansion, the first Brillouin zone (BZ) is converted 
into a region \x,y\ < 193/2 /Ltm [ll|. The population of 
the higher Bloch bands increases slightly during the lat- 
tice ramping down, but remains negligible. In fact, as 
shown in Fig.[TJc) and (d), there are few signals out of 
the first BZ. Although the initial states are quite differ- 
ent from each other, we find no distinctive differences in 
between the TOF images. We can understand this fea- 
ture by noting the following points. Since a wavevector 
k is not a good quantum number owing to the trapping 
potential, the difference between the momentum distri- 
butions of the AF and PM states is rather subtle even 
in the initial states. In addition, the effects of the lat- 
tice ramping down make the TOF images featureless, as 
discussed later. 

Here we discuss how to distinguish the AF state from 
the PM state using TOF images. For this purpose, we 
propose the following simple procedure: (i) before start- 
ing the lattice ramping down, we suddenly change the 
scattering length asc to asc = RscO'SC using a Fesh- 
bach resonance; (ii) then, start the lattice ramping down 
procedure as usual; (iii) at t(< iramp), before the proce- 
dure is completed, turn the trapping and lattice poten- 
tials off simultaneously; (iv) capture a TOF image after 
ballistic expansion in ixoF- We depict this sequence in 
Fig-ma)- We repeat sequence (i)-(iv) for different t val- 
ues, and investigate dynamics of fermionic atoms in the 
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FIG. 3: (Color Online) (a) and (b) Fraction of atoms in the 
2nd RBZ, F — 72nd/(-^ist +/2nd), as a function of t. The inset 
shows the definition of the 1st and 2nd RBZs. 



momentuni space. In the following, we adopt the same 
PM and AF states as in Fig.[T] as the initial states, and 
set tramp = 2 ms, and ixoF = 8 ms. Although Rsc is 
varied up to 10, the energy scale of the increased interac- 
tion strength is much smaller than that of the band gap, 
U = RscU < A. 

In Fig.[2l^b) and (c), we show TOF images obtained at 
t = 0.08,0.24, and 0.4 ms with Rsc = 5 for the PM and 
AF states, respectively. Even though the population of 
the higher hands is very small, all of the TOF images 
in Fig.[2{b) and (c) expand over the 1st BZ because the 
lattice ramping down is incomplete. We find that similar 
cross-shaped images are obtained a.t t — 0.08 and 0.40 ms 
for the AF state and at t = 0.08, 0.24, and 0.40 for the 
PM state. Significantly, a distinctive image is obtained 
at i = 0.24 ms for the AF state. We find that such dis- 
tinctive images appear periodically for the AF state from 
the calculations of a larger t. This kind of characteristic 
behavior cannot be seen for the TOF images obtained for 
the PM state. 

To clarify this characteristic behavior of the TOF im- 
ages for the AF state, we introduce a quantity; the frac- 
tion of the atoms in the second reduced Brillouin zone 
(RBZ) F = /2nd/(/ist + /2nd), whcrc list (/2nd) IS a sum 
of the intensity in the 1st (2nd) RBZ of a TOF image. 
The definitions of the 1st and the 2nd RBZs are given in 
the inset of Fig.[3l^a). The quantity i^ as a function of t 
for the PM and AF states with fixed Rsc = 5 is shown in 
Fig. 131a). We can clearly see the oscillation of F for the 
AF state at t < 1 ms, while we find no oscillation for the 
PM state. At t > 1 ms, F gradually decreases for each 
state. In Fig.[3l[b), we show the influence of Rsc values 
on the oscillatory behavior of F. The period of the char- 
acteristic oscillation for the AF state becomes faster with 
increases in Rsc- The estimated values of the period are 
0.44 and 0.21ms for Rsc = 5 and 10, respectively. 

We can understand this characteristic dynamics of AF 
states as the signature of translational symmetry break- 
ing. In our scheme, the scattering length is instantly 
changed from its initial value via a Feshbach resonance. 
This induces some scattering processes as a result of 
the additional mechanical work imposed on the system. 
The dominant effect comes from the potential scatter- 



FIG. 4: (Color Online) (a) Time evolution of the quasi- 
momentum distribution of atoms with k = (0, 0) and (tt, tt) 
for the AF and PM states, (b) Kinetic energy -Ekin, trapping 
potential energy Spot, and interaction energy Eint as func- 
tions of t. Thin lines (thick lines with plots) are data for the 
PM (AF) state. 



tial scattering is rewritten as S ^ ^-^^/X^ko- "^ku'^k+Qo-, 
where Q = (tt, tt) is the ordering vector of the AF 
states, and M is an averaged value of the magnetiza- 
tion Mr — l/2|(nrf) — («ri)|- This scattering, which 
transfers a wavevector Q, induces excitations from the 
1st RBZ to the 2nd RBZ, resulting in the above men- 
tioned oscillation. In contrast, for paramagnetic states 
whose translational symmetry in not broken, S does not 
transfer any wavevector. 

To discuss this oscillation in more detail, we calculate 
the quasimomentum distribution nk just before the bal- 
listic expansion. In Fig.UJ^a), we show the time evolution 
of rik=(o.o) and '^k=(7r,7i-) for the PM and AF states with 



fixed Rsc = 5- We find that, at t < 1ms, nk= 



.(0,0) 



and 



ing5 = f/Er 



(rira). If the density profile has the 



checker-board structure shown in Fig.[TJb), this poten- 



"k=(7r.Tr) for the AF state exhibit an oscillation with the 
opposite phase, clearly suggesting an excitation trans- 
ferring Q = (7r,7r). This oscillation results from the 
difference between the energies of quasiparticles with 
k = (0,0) and {tt,tt). In the AF state, the quasipar- 
ticle dispersion is approximately given by ek(k+Q) = 

2 f eok + eok+Q ± \/(eok - eok+q)^ + 4(A/f/)2 j , where 
eok — —'^J^xCos{kx)- Thus, we obtain the period of 
this oscillation as iosc ^ hTT/^J{W/2Y + {MUf. We ob- 
tain M ~ 0.3 from this equation, which is consistent 
with the value directly calculated from {n-ra) shown in 
Fig.[ljb). This feature allows us to estimate experimen- 
tally the magnetization value via an oscillation period. 

In Fig. Ufa), we find that n]^=(ofi) for both states in- 
creases at i > 1 ms, corresponding to the decrease in F 
shown in Fig.[2Ua). To clarify this behavior, in Fig.Hl^b), 
we show the change in the internal energy during the 
lattice ramping down for the PM and AF states. We cal- 
culate the kinetic energy ii^kin — X^ka £ok(f^k(T), trapping 
potential energy Epot = Vt Xro- ''^("■ro-), and interaction 
energy E{at = U ^^^{ura) {rirs) ■ The increase in E'lnt in 
the vicinity of t = ms is due to the sudden increase in 
asc- We see a large kinetic energy gain for both states, 
resulting from the fact that J increases exponentially as 
Vq is linearly decreased [IJ . Note that this effect is more 
dominant for a larger initial Vq (we now set Vq = lAEr). 
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FIG. 5: (Color Online) Fraction of atoms in the 2nd RBZ F 
as a function of t for the PM and AF states. The increase 
rate of the scattering length is changed as Rsc = 5, 1, and 0. 



The large gain of -Bkin is accompanied by an increase 
in the number of atoms around the bottom of the first 
Bloch band [k ^ (0,0)], and a decrease in those around 
the top of the first band [k ^ (7r,7r)]. This makes the 
TOF images featureless as mentioned above. 

Here, we consider the special cases of Rsc = and 1, 
where the oscillation of F disappears even for the initial 
AF state. We can easily expect that, without a change 
in the interaction strength {Rsc = 1): the character- 
istic oscillation will never occur. In addition, when we 
tune asc within the noninteracting limit {Rsc = 0), the 
characteristic oscillation vanishes because the potential 
scattering S (x U is absent. In Fig. [31 we show F for the 
PM and AF states with Rsc = 1 and 0, and also that 
with Rsc = 5 for comparison. As shown in Fig.ISJa), the 
F curves for the PM state change little as -Rsc is varied. 
In contrast, from Fig.[5jb), we can see the distinctive F 
curves of the AF state. Significantly, we can find that 
F decreases faster for Rsc — 1 than for Rsc = 0. This 
results from the fact that the relaxation that occurs dur- 
ing the lattice ramping down is accelerated by the energy 
gap of the AF states, Aaf ^ 2UM. 

Finally, we discuss the behavior of the F curves for the 



paramagnetic Mott insulating state that we cannot deal 
with in our mean-field calculation. Since this state with- 
out translational symmetry breaking has an energy gap, 
we expect to find that F curves show no oscillation but 
show fast relaxation for Rsc — 1 and 5, while exhibiting 
no characteristics for Rsc — 0. 

In summary, we investigated a 2D optical lattice sys- 
tem consisting of '^''K atoms with two different hyper- 
fine states. On the basis of precise numerical simula- 
tions of the TOF imaging process with lattice ramping 
down, we investigated the nontrivial dynamics that oc- 
cur during this process. We clarified that the kinetic 
energy gain induced by the change in the Hamiltonian 
results in the relaxation from the top to the bottom of 
the lowest Bloch band. We also suggested that this re- 
laxation process would be enhanced in the gapped phase. 
Furthermore, we proposed a TOF imaging procedure for 
detecting AF states, where we suddenly change the scat- 
tering length using a Feshbach resonance before starting 
the lattice ramping down and then investigate the dy- 
namical change in the TOF images systematically. In 
the translational symmetry broken states, the sudden in- 
crease in interaction induces potential scattering, which 
transfers ordering vector. This results in the oscillation 
of the momentum distribution of atoms, which can be 
detected with the standard TOF imaging technique. We 
showed that the signatures of AF states can be detected 
by the present method. Note that this method could be 
applied to other translational symmetry broken phases. 
In our calculations, we applied the mean-field approxi- 
mation to the interaction terms, and discussed only the 
dominant scattering effects. Although this might mean 
that we neglect the higher order scattering process, it is 
beyond the scope of our current work. 

We thank Y. Takahashi, S. Uetake, S. Taie, S. Sugawa, 
S. Suga and Y. Tokura for valuable discussions. 
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